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Abstract 

For the anisotropic K£iLi $) 2 + vJ2iLi -theory with N = 2, 3 
we calculate the imaginary parts of the renormalization-group func- 
tions in the form of a series expansion in v, i. e., around the isotropic 
case. Dimensional regularization is used to evaluate the fluctuation 
determinants for the isotropic instanton near the space dimension 4. 
The vertex functions in the presence of instantons are renormalized 
with the help of a nonperturbative procedure introduced for the simple 
gcj) 4 -theory by McKane et al. 

1 Introduction 



More than twenty years ago, Brezin, Le Guillou and Zinn- Justin (BGZ) stud- 
ied the phase transition of a cubic anisotropic system by means of renormal- 
ization group equations |jj . Within a (4 — e)-expansion, they found that to 
lowest nontrivial order in e, the only stable fixed point for N < N c = 4 is the 
0(iV)-symmetric one, where N is the number of field components appearing 
in the cubic anisotropic model. They interpreted this as an indication that 
the anisotropy is irrelevant as long as N is smaller than four. For N > A, 
the isotropic fixed point destabilized and the trajectories crossed over to the 
cubic fixed point. 
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Recently, our knowledge of perturbation coefficients of the renormaliza- 
tion group functions of the anisotropic system was extended up to the five- 
loop level by Kleinert and Schulte-Frohlinde |2|. Since the perturbation ex- 
pansions are badly divergent, they do not directly yield improved estimates 
for the crossover value N c where the isotropic fixed point destabilizes in favor 
of the cubic one. An estimate using Pade approximants @ indicates N c to 
lie below 3, thus permitting real crystals to exhibit critical exponents of the 
cubic universality class. 

For a simple 4 -theory, the Pade approximation is known to be in- 
accurate. In fact, the most accurate renormalization group functions for 
that theory have been obtained by combining perturbation expansions with 
large-order estimates and using a resummation procedure based on Borel- 
transformations 

It is the purpose of this paper to derive the large-order behavior of the 
renormalization group functions for the anisotropic g0 4 -theory. In a forth- 
coming paper we will combine these results with the five-loop perturbation 
expansion of Kleinert and Schulte-Frohlinde to derive the precise value for 
the crossover value N c . 

For the simple g0 4 -theory, the large-order behavior of perturbation coef- 
ficients has been derived by Lipatov |TI| [11], BGZ |Ll| and others ||13||— [15 



in a number of papers. The generalization to the 0(N) -symmetric case was 



given in |fL2| . An equivalent method for calculating the large-order behavior 
is based on the observation that for a negative coupling constant Green func- 
tions possess an exponentially small imaginary part due to the fact that the 
ground state is unstable |L6|, [L7|]- The imaginary part is associated with the 
tunneling decay rate of the ground state. It determines directly the large- 
order behavior of the perturbation coefficients via a dispersion relation in the 
complex coupling constant plane. 

In the semiclassical limit, the imaginary part of all Green functions can be 
calculated with the help of classical solutions called instantons. For a mass- 
less g(j) A -theorj in d — 4 space dimensions, these instantons can be found 
analytically. The imaginary part is a consequence of a negative frequency 
mode in the spectrum of the fluctuation operator, whose determinant enters 
the one-loop correction to the instanton contribution. McKane, Wallace and 



de Alcantara Bonfim [jnj found a way to continue the results of the g0 4 - 
theory in d = 4 to a field theory in d = 4 — e dimensions. They proposed an 
extended dimensional regularization scheme for nonperturbative renormaliz- 
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ing the imaginary parts of vertex functions. 

In the present work we have to extend this scheme to the case of a 4 - 
theory with cubic anisotropy, where the energy functional has the following 
form: 



H(4>) = J d d x 



N 



-v 2 



+ 



u 



N 



■2\2 



V N 

^ 8=1 



For N = 2, the corresponding model in quantum mechanics was first stud- 
ied by Banks, Bender and Wu [I9L |20fl who used multidimensional WKB- 



techniques to derive the large-order behavior of the perturbation expansion 
for the ground state energy. In 1990 Janke |^TJ presented a more efficient cal- 
culation using a path integral approach. In the present work, this approach 
will be generalized to quantum field theory and extended by a careful dis- 
cussion of the region near the isotropic limit v — > 0. This is important, since 
the infrared-stable cubic fixed point is expected to appear very close to the 
0(iV)-symmetric one. In fact, it will be sufficient to give the quantum-field 
theoretical generalization of [21] in terms of an expansion about the isotropic 
case in powers of v. 

The paper is organized as follows. The method is developed by treat- 
ing first the case N = 2. In Section 2 we derive the Feynman rules for the 
power series expansion of all Green functions around the isotropic limit. In 
Section 3 we calculate the small-oscillation determinants for the transversal 
and longitudinal fluctuations. In Section 4 we use the extended renormal- 
ization scheme of fllS] to find the full (real and imaginary) vertex functions, 



and derive the renormalization constants to one loop. In Section 5 we calcu- 
late the imaginary parts of the renormalization-group functions and thus the 
large-order behavior of the perturbation coefficients. In Section 6, finally, we 
extend the results to the physically relevant case N = 3. 



2 Fluctuations around the isotropic instan- 
ton 

For positive coupling constants u and v, the system defined by the energy 
functional (p]) is stable and the Green functions are real. On the other hand, 
if both coupling constants are negative the system is unstable and the Green 
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functions acquire an imaginary part. The corresponding functional integrals 
can be calculated by an analytical continuation from positive to negative 
coupling constants, keeping the factor exp [/(— j (Y^iLi 0;) 2 — \ Y^iLi 4>t) d A x] 
real. We perform this analytical continuation by means of a joint rotation in 
the two complex planes, substituting u — > u exp {i6) and v — > v exp {i6) , and 
rotating the azimuthal angle 9 from to n. At the same time, we rotate the 
contour of integration in the field space. The convergence of the functional 
integrals is maintained by the field rotation — > 0exp (— i8/A). 

A natural parameter for the anisotropy of the system is the ratio 5 = 
v/(u + v). The isotropic limit corresponds to 8 = 0. During the joint rotation 
of u and v, the parameter 5 remains constant. Thus, 5 is a good parameter 
for the anisotropy at both positive and negative couplings u and v. We shall 
use the coupling constants g = u + v and 5 for a calculation of the Green 
functions. These are given by the functional integrals 



G {2M \ 



Xi, X2-, 



'1,'2,->'2M 



JDcPexpi-H^]} 



(2) 



where the subscripts ik run through the iV components of the field <fi. As 
explained in the introduction, we shall first study the case N = 2 with the 
free energy functional 



>1,9>2 



d a x - 



1 r 



L (-v 2 
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:-v 2 



# + 2(1 + $ .(3) 



When expressed in terms of the coupling constants g and S, the Green 
functions possess an imaginary part for g < 0. For the reasons explained 
above, we shall derive an expansion of the imaginary parts of the Green 
functions (0) around the isotropic case, i. e., in powers of 5: 



oo A 

ImG=5> n W — 



p(n) 



exp (±\ [1 + 0{g)\ . 



(4) 
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Each power 5 n has its own n-dependent imaginary part. Given such an 
expansion, the large-order estimates for the coefficients of the powers g k 
follows from a dispersion relation in g: 

o 



G 



1 



7T 



,_1mG(g + iO) 
dg 

9-9 



(5) 



For a general discussion of the relationship between imaginary parts and 
large-order behavior see, for example, ch. 17 of [22 1. If the power-series 
expansion of G is 

oo 

G= E G kn g k 5 n , 

k,n=0 

the coefficients Gkn have the asymptotic behavior 

.A, 



G 



kn 



71 



■1) 



jfe!jfcP(»)-l [1 + 0(1/Jfe)]. 



(6) 



(7) 



These estimates apply to perturbation coefficients, in which the maximal 
power of g is much larger than the power of 5. Being interested in the 
region close to the isotropic limit, this restricted large-order estimation will 
be sufficient. 

An expansion of the Green function (0) around the instanton yields ex- 
ponentially small imaginary parts in both numerator and denominator. In 
order to isolate the imaginary part of the numerator, we simplify the denom- 
inator in (H) and calculate first the imaginary part of the approximate Green 
functions 



£(2M) ( 



Xl,X 2 , 



%2M, 



fD(f) ^(sQ&afo) • • • • • (j)i 2M (x2M)exp{-H[<f)]} 

JDcf>exp{-H o m 

(8) 

where the denominator contains only the free energy functional Ho. 

With the aim of calculating (|]) , we expand the fields around the classical 
solution of the isotropic limit 5 = for the space dimension d = 4. Thus we 
write: 

$= U L ^ C + U L i +U TV = ( ) 0e + ( C s ^^ ) e + ( ) V ' (9) 
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where <t> r . is the well-known qf(/> 4 -instanton in four dimensions 

(10) 



1/2 A 



\—g) 1 + A 2 (x — xq) 

and d is the rotation angle of the isotropic instanton in the (0i, 02)-plane. 
The fields £ and rj correspond to the degrees of freedom orthogonal to the 
rotation of that instanton. The parameters xq and A are position and scale 
size of the instanton, respectively. Inserting the expansion (|§D in if(0i,(/> 2 ), 
we obtain the expression: 

H{fa, <h) = H{<f> le , 2c ) + H 1 + H 2 + H 3 + H,, (11) 

where Hi collects all terms in £ and rj of order i. They are given by: 
(1) Linear terms: 



g y/2 J [l + A 2 (x-x ) 2 ] 2 
!) sin 2 (2tf) 8(2)^ / d d x- — — £ 



{-g)V* v ' w J [l + A 2 (x-x ) 2 ] 3 
+ _A_ smW 4(2)^/ A |i + A2( ^^ )2]3 ,, (12 ) 

(2) Quadratic terms: 

Hi((,v) = \f d d xZM L t+ l -f d d xr,M T r, 

+ 6^ m6 J d * X [1 + x2 £_ Xo)2]2 tr,, (13) 

where and are the operators 

^2 24A 2 „, „ 2 8(1 -5)A 2 
Mr = -V , M T = -V 2 ( 14) 

L [l + A 2 (a;-a; ) 2 ] 2 ' T [1 + A 2 (a; - x ) 2 ] 2 1 1 
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(3) Cubic terms: 

Hs(Z,l) = -(-8g)^jd d x 
sin 2 (2$) 



A 



i + x 2 (x — x y 



-5 
-5 

(4) Quartic terms: 



2 

sin(4$) 
4 



(-8^) 1/2 d 



A 



x- 



l + A 2 (x-x ) 2 
1 + A z {x — Xq) z 



(16) 



The linear terms ( 12[) contain factors e or 5, since the expansion of the 
fields around the isotropic instanton is extremal only in the four- dimensional 
isotropic limit. 

The classical contribution of the instanton is 



H( 



>lc, (P2c 



(17) 



A £ 8vr 2 



-(2 + lnvr + 7 ) 



X £ 5 87r 2 sin 2 (2^) | Q (5^ 



y 



The first term in i7(0 lc , 2c ) is evaluated up to the first order in e, because 
the one-loop renormalization will require replacing 1/g — > l/g T + 0[f(6 T )/e), 
where g r is the renormalized version of the coupling constant g. A contribu- 
tion of order e is not needed in the second, (^-dependent term, where it would 
produce a further factor 5, and thus be part of the neglected terms O(g) in 

(§■ 

Due to the anisotropy of the action, the fluctuation expansion (|12"|)-(|T7|) 
contains independent parts. However, all these terms are of order 5 and can 
therefore be handled by straightforward perturbation theory. 

Note that the angle •d disappears when the isotropic instanton is directed 
along the coordinate axis. Then the isotropic instanton coincides with the 
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exact solution in d = 4 for 5 > 0. This is also seen by inspecting the potential 
in (H). For 5 > 0, the term <pf + (p\ is larger than 2(1 — «5)0f 0|, so that the 
"tunneling-paths" of extremal action are obviously straight lines along the 
coordinate axis. 

For the region 5 < 0, the exact instanton follows from the known fact 
that the action (0) is invariant under the orthogonal transformation: 



01 = (01 + 2 )/V2, 2 = (0!-0 2 )/V2, 

with the new coupling constants: 



g=(2-6)g/2, 5 



25 



2-5 



satisfying 5 < for 5 > 0. 

In contrast to the method in [H], the treatment of the fluctuations in a 



power series in 5 does not allow us to deal with all modes perturbatively. 
Near the 0(2)-symmetric case, the Gaussian approximation for the rotation 
of the instanton becomes invalid, and the Gaussian integral must be replaced 
by an exact angle integration. The separation of the rotation mode must be 
done with the help of collective coordinates. The Jacobian of the relevant 
transformation can be deduced from the isotropic system. The field r\ of small 
oscillations must not contain modes of with eigenvalues of the order 5, 
since these would vanish for 5 = 0. The discussion of the longitudinal part 
is given in [^3[ . 

In dimensional regularization, all eigenvalues of which would be zero 
for d = 4 are of order e in 4 — e dimensions. To avoid eigenvalues of the 
order 1/e in the propagator for the longitudinal fluctuations, all collective 
coordinates of the four- dimensional case must be retained in d dimensions. 
Therefore the field £ in @ contains no modes of with eigenvalues pro- 
portional to e. 

In order to calculate the fluctuation factor and to separate the almost-zero 
modes from detMx, and detMr, it is convenient to do a conformal transfor- 
mation onto a sphere in d + 1 dimensions leading to a discrete spectrum for 



the transformed differential operators and fl24j , Pq| . Their products 
of eigenvalues can be given in terms of the Riemann (^-function which we 
easily expanded near e = 0. Simple 1/e-poles, which are characteristic of di- 
mensional regularization, appear directly from the known singularity of that 
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function. In place of the fields £ and 77, we define the fields and $2(p) 

on the unit sphere in d + 1 dimensions: 

$^0-!-^ $ 2=a l^/2 r/; ( lg ) 

where cr = 2/(1 + x 2 ). The instanton is supposed to be centered at the origin 
and rescaled by A. The spherical operator corresponding to the differential 
operator V 2 is 

V = \L 2 - \d{d - 2) , (19) 

where L 2 = J2 a ,b(.Padb — Pbd a ) 2 is the total angular momentum operator on 
the (d + l)-dimensional sphere. 

The eigenfunctions of Vq are the spherical harmonics Yj^(p) in d + 1 di- 
mensions They satisfy 

Vo Y l m {p) = -(/ + \d - + id) F^(p) , (20) 
and have the degeneracy 

(21 + d- l)T(l + d- 1) , s 

After the transformation, the angle-independent quadratic part of H 2 in (|T^) 
becomes 

l - J dm^-Vo - 6)$! + ~ J dfi$ 2 [-y - 2(1 - 5)} $ 2 , (22) 

where dVt is the surface element of the unit sphere in d + 1 dimensions. 

After rewriting the entire energy functional (|TT|) in terms of the new fields, 
we can summarize the Feynman rules for the diagrammatic evaluation of the 
functional integrals as follows: 

(1) Propagators (from the -^-independent part of if 2 ): 
longitudinal propagator: 

= "(-H-6)'] -1 (23) 

transversal propagator: 

-i 



{[-y - 2(1 -*)]'}" (24) 
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The prime indicates that the almost-zero modes are omitted when forming 

the inverse. 

(2) Tadpoles (H^: 



.£ ^AV /2 (T -l+ £ /2 $1 



1/2 



5 sin 2 (2^) ( r £ / 2 $ 1 



"2 J 



1/2 



(5 sin(4tf)<7 £/2 $ 2 



The bold dot stands for the ^-dependence of the vertex. 
(3) Two-point vertices (^-dependent part of H 2 ): 



(25) 



+ 



;sin 2 (2^)5($?-$ 2 ) 



sin(4#) S $i$ 2 



(4) Cubic vertices (H 3 ): 



(26) 



-{-2gf/ 2 {l-5)a- £ l 2 ^ 



+ 



A 



-(-2gf 2 5 Sin2(2 ^ } a- £ / 2 (3$ 2 $! - 

_(_ 2 ^) V2 5 S j^l a -e/2 ^3 _ 3$2$ 2 
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(27) 



(5) Quartic vertices (H4) 




+ 



+ 



X + X 



I 5 sin(4tf)cr 



$!$2 - $2^1 



x + x + x 



" sin 2 (2tf) £ 



-Z-6 
4 



(28) 



In order to obtain the leading corrections we must consider all connected 
diagrams of order 0(1/ g), 0(8/g), and 0(1). The contributions of these 
diagrams are added to H c in (|17D . We consider first the contributions to 
0(1/ g), which arise from the connected tree diagrams 



x — X, 




(29) 



and their $! insertions of zeroth order in g: 



X 



X X 



XXX 



(30) 



Since each <3>i-vertex produces an e-factor, the smallest order in e is two. 
However, terms of order e 2 /g are negligible for our calculation even after 
one-loop renormalization, since: 1/g l/g T + 0[f(5 T )/e] and e — > 0. The 
diagrams of 0(S/g) can be generated from those of 0(l/g) by one of the 
substitutions: 



X- 



(31) 
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By inspection, we see that the only diagram with an e-power smaller than 
two is given by 

X • (32) 

which is of order e 5/g. After the renormalization, a term of order 5 r appears. 
But this term can be neglected in comparison with the <5 r /g r -term from H c 
in expression (|ITD . Hence all the tree diagrams would enter only in a higher- 
order calculation. 

Contributions to order g° can appear from one-loop diagrams. The only 
possible candidates are 

^ + x — {' \ 

(33) 

where the 1 / e-pole from the loop integration and the ^-factor of the $i-vertex 
cancel. Following the method of [2~3|| , we have found that these diagrams 




contribute to the coefficient of the imaginary part of G^ 2M ^ in (H) a factor 

c{c\ = exp(-5 + 5 - 5 2 /2) . (34) 



3 Quadratic Fluctuation Determinants 

The angle-independent quadratic form is decomposed into a longitudinal and 
a transverse part. For the longitudinal fluctuations we obtain the determi- 
nant 

II 

<*•' 1=0 



' det V L ' 



det(-V r - 6) 
detf-V? 



{l + \d-d,){l + \d + 2) 
[l + \d - + \d) 



(35) 

which coincides with that appearing in the one-component g4> 4 theory, and 
is therefore known |2"3fl . It contains a bound state at / = 0, this being 
responsible for the expected imaginary part of G^ 2M \ and d + 1 almost-zero 
eigenmodes of order e associated with the dilatation and translation degrees 
of freedom of the instanton in 4 — e dimensions. Extracting these modes 
from the product 



f5D in the framework of collective coordinates |27|, 
obtain the well-known formal replacement rule for the determinant 



' det V L ' 
det Vol. 



J Vl (2 



7T 



we 



(36) 
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with 



3 3 7 
_ _| 7 

e 4 2 1 



iC(2) 



7T 



4 = (27r)- 1 / 2 5 5 / 2 exp 
and the Jacobian of the collective coordinates transformation: 

16vr 2 A £S 



jV L = X ( d -i) 



d+l 

2 



15g 



[l + 0(e,g)\. 



(37) 



(38) 



The expression for contains the Euler constant 7, and the derivative of 
the Riemann-zeta function ((x) at the value x — 2. Note the simple pole in 
e, which is related to the ultraviolet divergence on the one-loop level. 

The transverse fluctuations contain neither negative nor zero modes for 
5 larger than zero. The corresponding fluctuation determinant is given by 



' det V T x 
det V 0T , 



det [-V - 2(1-5)] 



detf-Vn 



n 

/=0 



(I + ±d - + ±d) - 2(1 - 6) 
(l + \d-l)(l + \d) 



■ (39) 



Just one zero eigenmode appears for I = in the isotropic limit <5 — 0, due 
to the rotational invariance in that case. For 5 > 0, the numerator in ( p9|) 
contributes a factor 1 / y/26. To avoid this artificial zero-mode divergence for 
5 — > we separate, as announced earlier, the angular degree of freedom from 
the integral measure via the collective-coordinates method. Similar to the 
longitudinal case, we make the formal substitution 



with 



' det V T ' 
detVbr, 



J Vt (2 



7T 



2 Co 



167r 2 A £ ' 
3^ , 



(40) 



(41) 



which coincides with the Jacobian of the isotropic model. Since the integra- 
tion interval for the angle $ is compact, no singularity appears in the limit 
5^0. 
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It is useful to illustrate the appearance of a divergent factor l/v25 in 
a careless use of Gaussian integral. Assuming 5 > 0, we expand the angle- 
dependent classical action around i? = up to quadratic order. The ensuing 
Gaussian integral can be evaluated after replacing the finite integration in- 
terval t? G [— 7r, 7r] by the noncompact one d G [—00,00]. In this way, the 
integral 

3vr 2 5 sin 2 (20)" 



\9\ 



is evaluated to 




eZ^exp 



1 /16tt 2 A £ 



3bl 



25 U 2 



25 



showing the spurious would-be zero-mode singularity, as the consequence of 
false rotation-mode treatment. 

Excluding the / = -mode in the numerator on the right-hand side of 



9|), we obtain for qf: 



2 l/2 



exp 



1 ~ (2Z + d - l)r(Z + d - 1) 
• > ■ - - 



{=i 



r(d)r(z + i) 



1 ~ (2Z + rf-i)r(f + d-i) 

2 f-r r(d)r(z + i) n 



1 + 



(Z + f + !)(/ + |- 2) 
(Z + f-l)(Z + f) 

25 



(Z + f + !)(/ + § -2) 



(42) 



The first sum is the well-known contribution from the isotropic limit. The 
second sum is to be expanded in powers of 5. For large Z, the sum diverges 
as e tends to zero. In order to separate off the divergence as a simple 1/e- 
pole, we use the zeta-function regularization method described in f23fl . A 
somewhat tedious but straightforward calculation leads to 



,-,T 



(2tt)- 1 / 6 3 1 / 2 



cxp 



(43) 

The series in front of the ^-singularity ends after the quadratic power of 5, 
thereby leading to the exact (1 — 5) 2 /3-coefficient of the simple e-pole. 
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Collecting all contributions to the one-loop expression of \m.& 2M \ we 
obtain the imaginary parts 

Im G {2M) (x u x 2 , • • • , x 2M )i 1 ,i 2 ,...,i 2M = 
-J d\d d x ^c^c L 2 c^2n)-^(2n)- 1 / 2 J v -J v -ex V 

x Il^(^) o / d$u Lh ($)uLiM ■ ■ -ULi 2M {$) exp[asin 2 (2tf)] \ (44) 



with 



A £ 8vr 2 S 



A = 



8tt 2 



1 - -(2 + hi7r + 7) 



(45) 



3g 2' 3 

and %k = 1,2. The expression contains a standard factor 1/2, since for 
symmetry reasons each saddle point contributes only one half of the Gaussian 
integral. 

It remains to perform the collective coordinates integration over the di- 
latation (A), translation (xo), and the rotation (■#) degree of freedom. Having 
obtained the imaginary parts of the Green functions, we go over to the bare 
imaginary part of the vertex functions. Taking the Fourier transform and 
excluding the (27r) d S(J2i (fo) -factors, an amputation of the external legs of 
G { - 2M \x 1 ,x 2 , . . leads to: 



ImT 



(2M) 



(?m) 



H,l2,---,l2M 



x \ J d#u Lil ($)u Li M ■ ■■u Li2M (#) exp[asin 2 (2tf)] [1 + 0(e,g)} (46) 
o J 



with 



0(g)=2 d / 2 7 r^3 1 / 2 |g| 1 - (rf/2) K f „ 1 (|g|), (47) 
where K^(z) is the modified Bessel-function. The coefficient Cb is given by 



C b = 2- 2 / 3 3 1 / 2 7T- 11 /3 exp 



l(l0-25 + 5 2 )+^P-4 7 -4+^ 
3e v ' n 2 9 
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(48) 



An e-singularity results from the ultraviolet divergence of one-loop diagrams. 
As discussed in the following section, this simple pole in e is removed by a 
conventional coupling constant renormalization. 

The contribution from the imaginary part of the denominator in equation 
(H) follows for M = in fl46|) . We observe that it is of higher order in g in 
comparison with the imaginary part of the numerator, and can therefore be 
neglected. So the imaginary part of (H) supplies the desired imaginary part 
of®. 

4 Renormalization 

In the preceding section we have calculated the contribution of the quadratic 
fluctuations around the isotropic instanton. The almost-zero eigenvalues of 
translation, dilatation and rotation have been extracted and treated sepa- 
rately. Of course, the resulting expressions are useless, if we are not able to 
renormalize the theory. This is the most difficult additional problem which 
arises in the transition from quantum mechanics to higher-dimensional field 
theories. 

A systematic scheme to renormalize both the real and imaginary part 
of vertex functions for a g0 4 theory in 4 — e dimensions was introduced by 
McKane, Wallace and de Alcantara Bonfim. They calculated the full (real 
and imaginary) renormalization-group constants using an extended minimal 
subtraction scheme to one loop (the conventional one is given in |5D|). 
We have extended their method to the case of two coupling constants in view 
of applications to critical phenomena. 

For an illustration we consider first the four-point vertex function. The 
result fl46|) of the one-loop calculation about the instanton reads: 

ImT^ 4 (q m )i,j,k,i — 
00 f (S/2) n 1 2 r 

£ -^T^o / d^uU^n L ^)u Lk {d)uU^Hm] 2n 

ri I 77/. — J 
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-c b 



dX 

T' 



A £ A 



rf + li + 2». 



exp 



A e A 



n 

i/=i 



A , 



(49) 



with A = + 0(e) and z, j, k,l = 1, 2. It remains to evaluate the integral 
over the parameter A. This is a relict of the introduction of collective coor- 
dinates in the instanton calculation. The integral converges for small A due 
to the exponentially decreasing of modified Bessel-function Kd/2-i- For large 
A, the product (^/A) 2 0(g l/ /A) behaves like 



2 2 3 l/2 



7T < 



1 2A 



<h_ 
2A 



l + -e (7-Iutt) + C( £ 2 ) 



+ l-7-lnvr + C>(e; 



*1 
2A 



(50) 



so that the integral diverges logarithmically as e goes to zero. This divergence 
causes a 1/e-pole in the imaginary part. Since after the above approximation 
of 0(g„/A) the integral diverges for A — > 0, a small A cutoff \i has to be 
introduced. In this way, the A-integral in ( fi9|) takes the form 



Cb 



dX 



X s A 
9 , 



d+6+2n 



exp 



X £ A 
. 9 , 



(2 2 3 l/ 27r) 4 ^ ^ ^ 



e>( — 

2A 



(51) 

The coupling constant g is chosen to lie on top of the tip of the left-hand cut 
in the complex g-plane. Being integrated only near g — > 0~, the integral can 
be evaluated perturbatively in g. Using 



J^A a£ A £ exp -( 

it L \ 



A e v4 



1 M or 

— ^V ae exp 



[1 + 0(g)] , (52) 



we get for ( pT] ) the expression 



8q2 4^ 1 # 
Z O 7T Cb — — 

£ A 



d+6+2n 

H £ A\ 2 



exp 



S 5 



(53) 
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After expanding the factor (q v /n) e in powers of s, we obtain the typical 
finite contribution \n{q v /fi). Since this term can be omitted by a minimal 
subtraction, we are left with 



1 oo 
£ n=0 



-2sy 



n!r(2n + 3) \A 



d+6+2n 

H e A\ 2 (if A" 

— exp It. 



where 



and 



(54) 



1 , i = j = k = I 
, otherwise . 



(z,j,M = l,2) 

The tensor structure arises explicitly upon doing the integrals 

T(a + l)r(/3 + 1) 1 



tt/2 

i sin 2Q+1 (tf) cos 2/m (tf) dtf 



2 T(a + /3 + 2) 



B(a + 1,0 + 1) 



for all combinations of indices in Im T^].;. The first term in ( |54"D contains only 
mixed index combinations, whereas the second term has only a contribution 
from equal indices. 

Now we can proceed to renormalize the bare results. The real part of 
T\f kl (q) is the perturbative one, and it is easily calculated to one loop. In 



terms of our coupling constants, the real parts ReT-^(g) take the form 



Rer[ 4] \q m )i,j, 



k,i 



-6(1-% + 
3 



+ -60 



4ir 2 e 



ijk 



^ (10 - 145 + AS 2 ) + O(e )^ V £ } {S ijkl - 5 ijkl } 
10 - 25 + 5 2 ) + O(s ) g 2 ^ S m , (55) 
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where /x is the arbitrary momentum scale introduced above in (pi]) . 

In the absence of an instanton, the wave function renormalization has no 
one-loop contribution. In the consequence, the expression Q55D is rendered 
finite by a coupling constant renormalization only. The subscript r of the 
coupling constants indicates renormalized expressions in the conventional 
perturbative manner, i. e., 



9r = 



- -V (10 - 28 + 5 2 ) g y 



2e 



5 V = 5 + 



]ir 2 e 



-25 + 5 2 + 5 3 ) gfi- £ 



(56) 



Inserting this into (p4|) and ( |55l) we find the perturbatively renormalized 
vertex functions 



Fi 4 \q m )i t j,k,i = -6/i £ [Fi(5 T , g T ) {S ijk i - 5 ijk i) + F 2 {5 r , g v )5 ijM ] (57) 



with 



Fi(5 T ,g T ) = 



£ n f„ n! T(2r! + 3) 



) d+4+2n 
exp — 
\9r, 



(58) 



and 



^2 (S T , gr) = 

1 00 (-2S ) n r fn + |) r fn + | 
S r +*2 7 37r 4 c r -£ 1 rj 1 2/ 1 " 



T(2n + 3) 



ii+4+2n 

— exp — 

9r XQr, 



(59) 



The coefficient c r is given by 

A 



c r = c b exp 



10 - 2<5 r + 5; 
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3l/2 2 -2/ 37r -2 exp 



2<5 r + (5; - — (10-25 r + 5, 



3s 



4C'(2) 7 2 



3 7 



- + 0(5 r ) + C»(5 r ^ 



(60) 



The 1/e-pole terms cancel. Thus, c r remains finite for e — > 0. For the 
leading expansion (f|) we can take the perturbative renormalized expres- 
sion ( |60D at the position 5 T = 0. The only remaining singularity is the 
1/e-factor in Imr^^g)^. It requires a further renormalization. In our 
special choice of coupling constants, the disastrous fc!-divergence of pertur- 
bation series appears in the expansions coefficients of g k . Since 5 is a di- 
mensionless anisotropy measure for positive as well as negative g, we have 
expanded the imaginary part of vertex functions around the isotropic case in 
a well defined manner. For the derivation of nonperturbative corrections to 
the renormalization-group functions and for later applications, however, it is 
convenient to avoid ratios of coupling constants. Therefore we proceed with 
v r = g r <5 r and g r instead of 5 r and g r . Then v T receive the rule of <5 r . 

Now we apply the extended minimal subtraction and perform a second 
(nonperturbative) renormalization step to eliminate the 1/e-pole term: 



9k = g t + i2 7 3n A c t -^ i 



n=0 



< / 3 
n! V47T 2 




n n ( 1 5 
B [n H — , n H — 



exp 




Y oo 

v r — i2 7 37r 4 c r -y^ 

£ n=0 L 




(61) 



The subscript R denotes the fully renormalized coupling constants. The 
required wave function renormalization in nonperturbative terms give no 
leading contributions to the imaginary part (see below). Therefore the 1/e- 
singularity in (|57|) is removed by a nonperturbative coupling-constant renor- 
malization only. 
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We now proceed to investigate the two-step renormalization for the wave 
function. There are two additional complications for the -functions. First, 
T m r( 2 ) is quadratically divergent. Second, there is a momentum- dependent 
divergence in the imaginary part of T^ 2 \ For the g0 4 -theory it was shown 
in that the undesirable momentum-dependence disappears during the 
nonperturbative renormalization process. This is comparable to the situation 
in conventional two-loop perturbation expansion. We continue to show that 
this non-obvious cancellation still works in the more complex case of two 
coupling constants and more than one field components. 

According to ( [46] ) the imaginary part of the unrenormalized two-point 
vertex function reads 



Imr( 2) (g),. 



' 1 T 



A 



X s A 



d+4 
2 



exp 



2rr 



X- / d$u Li ($)u Lj ($)exp asin 2 (2$) }. (62) 



The integration over the angle $ can be done: 



2~ 



~ / d$u Li ($)u Lj ($)exp asin 2 (2$) 



r(n + |)r(n + 



n=0 



nl 



T(2n + 2) 



n=0 



nl 



1 3 
n H — , n H — 
2' 2 



(2) 



(63) 

Thus, we get the usual form: Im I^ 2 - • = 5ij Im . The wave function renor- 
malization is contained in -£pT^ 2 \q). Following a similar procedure as in the 
case of r^ 4 ), we obtain the imaginary part 



2V3c b -£ 



-2sy 



n=0 



nl 



r>, 1 3 

B \n + -,n + - 



x 



H £ A' 



d+2+2n 
2 



exp 



(64) 
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Apart from the appearance of the (1/e) ln(g//i) term, a very similar expres- 
sion is found for ImT^\q). In order to show the cancellation of this term 
during the second stage of renormalization we need the two-loop expression 
for the real part of ■Jgz'T™ (q) 



4(8tt 2 ) 2 £ 



+ P + 0(e) 



(65) 



where P denotes some number which is the contribution of order e°. At 
q = fi, (|65|) is commonly defined as the inverse wave function renormaliza- 
tion constant (Z^)^ 1 . In terms of the perturbative renormalized coupling 
constants is given by 



Z P 



1 - 



(4 - 26 r + 5 r 2 ) g 2 



and together with 
d 



4(8vr 2 )% 

we find after the first step of renormalization 



(66) 



— T?\q) = 1 + ( ~9r - 77^r + -V] 



4 



1 00 r n n 
+ *2 4 vr 2 3c r -y 



3 

< / 3 



2(8tt 



2^2 



V47T 2 



1 3 

B \n + -,n+ - 



A 

9r, 



d+2+An 



A 

exp — 



(67) 



We have again introduced the coupling constants v r and g T which are more 
convenient for practical applications. This expression reads in the fully renor- 
malized couplings (|6~ID 



_d_ 

dq 2 



r«(?) = i+(^-|v a «R + ^) (p 



2(8n 



2\2 



In 



+ z 2 4 vr 2 3 c r -J2 



n=0 



111 



4vT 2 



X 



2 4 tt 2 A 5 ( n + -,n + 



5\ 4(n + l) 
2n + 3 



2(8vr 2 ) 2 



In 



q 
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X 



+4 i +5-« + DM- hi2+in u 

d+2+4ra 

9r) ^ \9r) J ' 



(68) 



By using the definition of the I?-function in terms of Gamma-functions one 
can easily read off the cancellation of (1/e) ln(g///) singularities for every 
power n of the coupling constant Therefore we are left with 



d 



— lf> (?) = 1 + ( -<? R - -wr(?r + - 



4 



R 



2 8tt 



2\2 



In i 



+ z 2 4 7r 2 3c r ^ (Vtt 2 A P + 7- - - In 2 



oo 
n=0 L 



fS / 3 \ n / 1 3 

-4 — o B (n + -.n+ - 
n UttV V 2' 2 



) d+2+4n 
exp — 



(69) 



Now we can define the nonperturbative minimally subtracted wave function 
renormalization to remove the 1 / £-pole: 



Z 4> 

np 



l-i 2 4 7T 2 3c r ^ (^2 4 7r 2 AP + 7-i-ln2 



oo 

*E 

n=0 



V R f 3 \ r, f 1 3 

-r — n B (n + -,n + - 
n! UW V 2' 2 

d+2+4n 

x exp — 

V 9rJ \9rj 



(70) 



The full renormalization constant is defined by Z^ := Z£ and follows 
from the product of flBTf ) and (|70|). 

A further important Green function is associated with the composite field 
(l/2)0 2 : 

G {2 ' 1) (x 1 ,x 2 ;x 3 ) iJ = (^ i (x 1 )(j) j (x2)-(t)l{x 3 )^ . (71) 
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By a straightforward application of the techniques discussed in the preceding 
sections we obtain for the corresponding vertex function 



Imr < £' 1) (q 1 ,q 2 ;q 3 ] 



1 

Ck 

2 



dX 



d+6 



exp [~, 



9 



2 



a; 2 



2tt 



d$u L 0)u L j{$) exp a sin 2 (2$) 



where 



0(g) = 2^ 2 )/ 2 vr^)/ 2 3 |g| (4 - d)/2 ^ 4)/2 (|g| 



(72) 



(73) 



Now we can proceed as in the case of wave function renormalization. Per- 
forming the collective coordinates integration over A and •& we obtain 



Imr^' 1] (q 1 ,q 2 ;q 3 )i, j = 5y 2 4 vr 2 3 2 c b - ^ 



1 



'-2 6) n 



n=0 



n! 



x 



which agrees with the usual form T^ 1 - 1 = 5^ T^' 1 ' . At q 3 = 0, q\ = q\ = /i 2 , 
the vertex function T^ 2 ' 1 \q 1 ,q 2 ;q 3 ) defines a third renormalization constant, 
called Z^, via 

r ( 2 ,i) = (z+zfy 1 . (75) 
The perturbative Z^ 2 follows from the real part of I^ 2,1 ); 



d+4+2n 

H e A\ 2 



B[n+-,n + - 



(2,1) 



exp 



(74) 



(76) 



where Q is a number of order e°. In terms of the perturbative renormalized 
coupling constants, the one-loop expression of Z% is given by 



1 



Z t =1 + ^( 4 -^. 



8tt% 
24 



(77) 



Including Rel^ and Imr^ 2,1 ^ we obtain after the perturbative step of renor- 
malization 



3* - 



Q 



1 00 r v n 



8tt 2 

x ( 3 \ n 



x 



V47T 2 

' ^4 s 

, ft. 



1 3 

B [ n + -,n + - 



d+4+in 
2 



exp — 



(78) 



In the second step of renormalization the (1/e) ln(g//i) singularities cancel, 
and we find 



1 /8tt 2 



-«2V3'c r - — Q -7 + ln2 



00 

*E 

n=0 L 



£ \ 3 

3 ^ n 



n\ \Atc 2 



D [11 + -.11 + -) ( — 



d+4+4n 



/ A 
exp — 

\9R y 



Hence, the nonperturbative Z^ 2 is 

ZtXl = 1 + i 2^ 2 3 2 c r - ( - 7 + In 2 



£ V 3 



(79) 



00 

xE 

n=0 



4(1 

ra! V47T 2 



/-> ( „ + _. „ + - ) [ - — 

K 9R, 



d+4+4n 



— A' 
exp — 

\9r, 



With respect to the leading expansion @) this is equal to . 



(80) 



5 Renormalization group functions 

After having calculated the imaginary parts of the renormalized couplings 
and of the renormalization constants Z^ and Z^ , we are in the position to de- 
rive the nonperturbative renormalization group functions /3(<7r, Vr), 7(fi'R, Wr) 
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and 7<£2 (<7r, vr) in the usual way. First we derive the /3-functions from the 
expressions fl6T|) for g R (fi) and fR(/i), respectively, and obtain: 



d_ 



d d 7 4 1 ~ 9 / 3 



1 5 

B(n+-, n + - 



t4 

9r, 



d+4+4n 



9 d 7 4 1 ~ d 

fJL— V R = V-Q^Vr ~l2 37T Cr-}^^ 



— ( A 
exp — 

> 3 



x 



n=0 ^ 

An 



n\ \47T 2 



2n + 3 



9r, 



d+A+An 
2 



B [n + -,n 



A 

exp — 



(81) 



The derivatives are taken at fixed bare couplings g and v. These functions 
will be finite as e — > 0. Furthermore, they will have a left-hand cut in the 
complex g-plane. On top of the tip of this cut we find 



AXflR, vr) = (3%(g v , v r ) - i2 7 3iT A c r J2 



n=0 



n\ V47T 2 



A" 



d+6+4n 



r, , 1 5 

B [n + -,n + - 



'A' 



x -- 



9x 



AX<?R, v R ) = /3»( gi , V t ) + ^2 7 37T 4 C r ^ 

n=0 



exp — 



X 



4n 



n! \47T 2 

v 9r 



1 5 

B[n + -,n+- 



d+6+4n 



,2n + 3 

The real parts are the familiar perturbative /3-functions 

1 



exp — 

\9r, 



(82) 



Pv(9r,V r 



~£9r + 



-sv r + 



3 



-2g T v r + vf + 10g; 



(4\g T v T 



(83) 
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In the leading contributions the pair of coupling constants (g T , v T ) may be 
replaced by (g R ,v R ) in (|52[) . This can be proven using fl6ip . The result is 



P g (g R} v R ) = (3 p g {g R ,v R ) -i2 7 3vr 4 c r ^ 



n=0 



^R / 3 



±M n + -,n + - 
n! V47T 2 / V 2 2 



x -- 



A_ 

9R, 



rf + li + 4n 



A 

exp — 



Pv(9r,vr) = P%(g R ,v R ) +i2 7 3ir 4 c T J2 



V R ( 3 



X 



n=0 

An 
2n + 3 



ryf 1 5 

, , B \n H — , n H — 
n! V4vr 2 y V 2' 2 



9k, 



d+6+in 



exp 



'A' 



(84) 



Note that the /3-functions have a form expected within a minimal subtracted 
scheme: f3 g (g R ,v R ) = -eg R + /3*(g R , v R ) and f3 v (g R , v R ) = -ev R + (3*(g R , v R ), 
where the /3 4 denotes the /3-functions in four dimensions. The critical large- 
distance behavior of the correlation functions is defined by the anomalous 
dimension of the field 4>, which is defined by 

r\ r\ r\ 

7 (<7r, v R ) = //— In = n— In Z* + ^— In Zt p . (85) 
With (|66| ) for Z^ and ([TDD for we get the ^-independent expression 



2 8tt 



l — (Agl - 2v R g R + v R ) + i 2 4 vr 2 3 c r (VvrV 1 P + 7 - ^ - In 2 



oo 

><E 

n=0 



^R / 3 



n! \47T 2 



d+A+An 

»+§)(—) 2 -p(£ 



.(86) 



The divergence of the length scale at critical temperature is governed by the 
anomalous dimension of the composite field </> 2 : 



lp (<?R, Va) = In ^ 2 = In Zf - ^ In z£ . (87) 
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Using ( fT7[ ) and (|80|), we obtain 
1 



(4^R-^R)+^ 2 4 7r 2 3 2 c 1 



57T 



■Q- 7 + ln2 



oo 
n=0 



d+6+4n 



v% ( 3 



0R, 



(88) 

For the leading contributions to the imaginary part, it was possible to replace 
(g r , v T ) by (#r, vr). We remark that the expansions , (|5E|) and (|8lf ) have 
the correct isotropic limit for ur — > 0. 



6 Generalization to the case N = 3 

Generalizing ([J) to iV = 3 we write: 



u L (p c + u L £ + u Tl r] + u T2 x 



— sin d 



( cos 9 cost? 



-h < ) + I cos$ I // + 



/ sin # cos d 
= sin 6 1 sin $ 
\ cos # 

where 9 and $ are the rotation angles of the isotropic instanton. This change 
of variables yields the following expansion of the energy functional: 



cos 9 sin d | X > (89) 
\ - sin 9 



H(<f>) 



1 f d d xiM L i d d xr)M Tl r) + 1 / d d x X y h, \ 



2 

eA{2f' 2 



d d x 



A 3 



(-g) 1/2 

(-8g) 1/2 I d d x 



[1 + A 2 (x-o;o) 2 ] 2 
A 



1 + A 2 (x-x ) 2 
+ 6 F(g,&,9,Z,ri,x) 



(90) 



with Ml and = Mr 2 = My of equation ([14]). For brevity, we have 
written down explicitly only the terms responsible for the leading contribu- 
tions in the expansion (Hf). The remaining terms are denoted collectively by 
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S F(g, 9, £, 77, x)- The classical contribution to the energy functional is 
, s A E 8tt 2 r 



H((f>i 



c, H>2ci <P3c 



9 3 
X £ 5 8tt 2 



1 - - (2 + lnvr + 7) 



2g 3 



sin 4 9 sin 2 (2t?) + sin 2 (29) J + I -e J . (91) 
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With similar steps as for iV = 2 we obtain the bare imaginary parts of the 
vertex functions: 



(3) f dXj . d _ M ( d _ 2 ) 

*■ 



A £ 8vr 
^9 



2 \ (rf+3+2M)/2 



cxp 



„ 1 n v N 



27T 7T 

x-|| mi9d9d'&u Lh {^e)u Li2 {-&,e)---u Li2M {-&,e) 




x exp 



a sin 4 9 sin 2 (20) + a sin 2 (26) \ [l + 0(e,g)] 



(92) 



where 



c (3) = 2 -5/6 37r -13/3 exp 



1 

3e 



11 -A5 + 25 2 ) + 



5C'(2) 9^ 17 
4 



7T 



2 2 7 



(4 = 1,2,3), (93) 



with a, A as in (fl5|). After evaluating the integrals over the collective coor- 
dinates A, 9, and •& we obtain for the imaginary part of the four-point vertex 
function 

I U1 ^b\ ( lm)i,j,k,l = 



2 8 3 2 7r 4 ct 3) 



rf+7+2n 

(-2 5) n fg\( fi £ A\ 2 ffi £ A 



e^ {T(2n + 7/2) \A 



exp 
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« r(3 + 2n-p)r (p+\ 

X E ' 



+r (n - p + ^ r (n - p + D % fcZ I . 



3r ( n - p + - ) (Sijhi - S ijk i) 



(94) 



Similarly the bare imaginary parts of -^F^(q) and r( 2,1 )(g 1? q 2 ; (fe) are found 



to be 



_d_ 

dq 2 



4^2o J3) 



7 -- + ln 



2/i 



1 00 

c n=0 



-2 5)' 



d+3+2n 



(95) 



and 



Imff g 2 ; g 3 ) M = 2 4 tt 2 3 2 - £ 



4_2 Q 2 (3) 



-2S) n 



n=0 



I(n) 



x 



-HI, 



rf+5+2n 

' exp — , 



9 



(96) 



with 



7 W = £( p ^B(2 + 2n-p,p+ 1 -^B(n-p+ 1 -, n-p+|). (97) 

Now we apply the extended renormalization scheme. Using the coupling 
constants t> and g, we obtain for the nonperturbative renormalized couplings 



Y 00 

g R = g T + t2 7 3n 4 cl 3) -J2 

£ n=0 



1 00 
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(98) 
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with 



J v {n) = E 

p=0 



n 
P 



4(n — p) 
2{n-p) + 3 



B ( 3 + 2n — p, p 



r, I 1 5 

xB [n — p + -,n—p+- 



(99) 



and 



c (3) = 3 2-5/6^-5/2 



exp 



5C'(2) 



7T 



(100) 



For the nonperturbative renormalization constants and Z^ 2 we find 



Z* p = l-* 2 4 7r 2 3c( 3 )-( 2 4 7r 2 AP + 7 -^-ln2 
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(101) 



and 



ZtXp = 1 + i 2 4 vr 2 3 2 <f 1 - 7 + In 2^ 
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xE 
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n! V47T 2 



exp — 
V#R, 



(102) 



To calculate the nonperturbative renormalization group functions /3(g R , v r ), 
7(<?r, i>r) and 7^2 (g R , t> R ), we simply repeat the calculations of the case N = 
2. Using the expressions (|98|) for #r and t> R yields the /3-functions 
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/%r^r)+22 7 37t 4 c( 3) E 
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n\ V47T 2 



J v (n) [- — 
V 9r, 



d+7+An 



exp — 
\9R, 



(103) 



The expressions for the nonperturbative anomalous dimension of the field 
and the composite field 2 are given by 

l(9R, vr) = 7 P G?R, vr) + i 2 V3 cf (2 V3" 1 P + j - i - In 2 
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,(104) 



and 



7^(sh,vr) = 7j 2 (SR^R)+2 2 4 7r 2 3 2 c( 3) ^-Q- 7 + ln2 
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exp — 



(105) 



respectively. 



7 Discussion and Conclusions 



Our main results are the real and imaginary parts ( 103|) -(105) of the 



renormalization-group functions /3, 7 and 7^2 in terms of the nonpertur- 
bative renormalized couplings g& and vr. Via the dispersion relation fl5|), we 
obtain the large-order behavior (^). It is useful to check the expansion terms 
proportional to powers of l/\g\ accompanying the v n in the imaginary part 
of each renormalization-group function. If the power of l/\g\ is denoted by 
p(n), we note that p(n) is the same as for the corresponding vertex func- 
tion before the integration over the dilatation parameter A. This integration 
with the measure din A reduces p(n) by 1 [see ([52])]. However this effect is 
compensated by the differentiation with respect to the logarithm of the scale 



32 



parameter // in the definition of the renormalization group functions. The 
power p(n) in the vertex functions can simply be explained: For a 2M-point 
vertex function with k 2 -insertions there is first a power (l/\g\)( M+k ^ from 
the fields. This is multiplied by a factor (l/|(y(|)( d +( Ar ~ 1 )+ 1 )/ 2 from the d trans- 
lational, (N — 1) rotational, and one dilatational would-be zero modes. The 
classical contribution to the energy functional yields a factor (l/|(j|) 2 for each 
power of the anisotropic constant v. Hence, for the renormalization group 
function which follows at the one-loop level from the vertex function Y^ 2M ' k \ 
the value of p(n) is given by 

d + N + 2(M + k)+4n 
p{n) = . (106) 

Inserting the imaginary parts into a dispersion relation (|5|), we can esti- 
mate the large-order coefficients of the corresponding series expansions 

oo 

f(9,v)= E hn9 k v\ (107) 

fc,n=0 

where n <^ k. The function / stands for /3 g , j3 v , 7 and 7^2, respectively. The 
rules how to go from the imaginary parts of / to are known from fl7|). For 
instance, the large-order coefficients of the /5-functions for N = 3 are given 
by 

( 3 \ n JJn) 



n! 
3 ^ k 



x (-l) fc (JL^J k\ k 2n+( - d+ ^ /2 [1 + 0{l/k)} 



Pv,kn = -4 3) 2 7 3tt 3 



3 \ n J v (n) 



An 2 J n\ 



3 x k 



x (_!)*(_!_ J k\k 2n+ ^/ 2 [l + 0(l/k)} , (108) 

with J, c£ 3 ) from (|99|) and ( |100| ), respectively. 

The leading large- A; behavior can be checked by means of a simple com- 
binatorial analysis. We start with a series in terms of the standard couplings 



33 



u and v [see ([!])]. For illustration, we consider the contribution of the pure 
li-powers. It is known from a theory with only one coupling constant, that 
the coefficients have the large-order behavior 

f lU l — ^{-a) l T{l + b+l)u l (Z>1). (109) 

Changing over to our couplings {u = g — v), the right-hand side of ( |109| ) 
contributes 

7(-iy- fc (-«yr(/ + 6 + i) ( l k ^jg k v l -\ (no) 

with 

/ i \ r(/ + i) 

\k ) k\{l-k)\' 
Replacing I by I = n + k, the contribution of g k v n has the form 

e kn i \n, \k ^(fc + n + b + 1) T(k + n + 1) fc 

An^S n ^7« — rrV L g k v n . (in) 

Now, for large k, the T-functions can be approximated by 

T{k + 5 + l) — > k\ k 5 (1 + 0{l/k)) . (112) 
Thus in the region n k we obtain 

f kn ~{-a) k k\k 2n+ \ (113) 

in agreement with ( |108| ). 

As stated in the introduction, the stable cubic fixed point is expected to 
lie in the vicinity of the isotropic fixed point. Since v becomes very small in 
this region, reasonable results should be obtained by resumming the ^-series 
accompanying each power v 11 . This will be done in a forthcoming publication. 
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